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Abstract 

We report on some recent results regarding the dynamical behavior of a trapped Bose-Einstein 
condensate, in the limit of a large number of particles. These results were obtained in a joint 
work with L. Erdos and H.-T. Yau. 

1 Introduction 

In the last years, progress in the experimental techniques has made the study of dilute Bose gas near 
the ground state a hot topic in physics. For the first time, the existence of Bose-Einstein condensation 
for trapped gases at very low temperatures has been verified experimentally. The experiments were 
conducted observing the dynamics of Bose systems, trapped by strong magnetic field and cooled down 
at very low temperatures, when the confining traps are switched off. It seems therefore important 
to have a good theoretical description of the dynamics of the condensate. Already in 1961 Gross 
[7J |H] and Pitaevskii ^1] proposed to model the many body effects in a trapped dilute Bose gas by 
a nonlinear on-site self interaction of a complex order parameter (the condensate wave function ut). 
They derived the Gross-Pitaevskii equation 

id t u t = -Au t + 8ira \ut\ 2 u t (1) 

for the evolution of ut- Here ao is the scattering length of the pair interaction. A mathematically 
rigorous justification of this equation is still missing. The aim of this article is to report on recent 
partial results towards the derivation of (^Q) starting from the microscopic quantum dynamics in the 
limit of a large number of particles. Here we only expose the main ideas: for more details, and for 
all the proofs, we refer to [I]. 
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Also in the mathematical analysis of dilute bosonic systems some important progress has recently 
been made. In ^3], Lieb and Yngvason give a rigorous proof of a formula for the leading order 
contribution to the ground state energy of a dilute Bose gas (the correct upper bound for the energy 
was already obtained by Dyson in [2J, for the case of hard spheres). This important result inspired a 
lot of subsequent works establishing different properties of the ground state of the Bose system. In 
|12| , the authors give a proof of the asymptotic exactness of the Gross-Pitaevskii energy functional for 
the computation of the ground state energy of a trapped Bose gas. In , the complete condensation 
of the ground state of a trapped Bose gas is proven. For a review of recent results on dilute Bose 
systems we refer to All these works investigate the properties of the ground state of the system. 
Here, on the other hand, we are interested in the dynamical behavior. 

Next, we want to describe our main result in some details. To this end, we need to introduce 
some notation. From now on we consider a system of A bosons trapped in a box A C R 3 with 
volume one and we impose periodic boundary conditions. In order to describe the interaction among 
the bosons, we choose a positive, smooth, spherical symmetric potential V(x) with compact support 
and with scattering length ao- 

Let us briefly recall the definition of the scattering length ao of the potential V(x). To define ao 
we consider the radial symmetric solution f(x) of the zero energy one-particle Schrodinger equation 

(-A + ±V(x))f(x) = 0, (2) 

with the condition f(x) — > 1 for \x\ —* oo. Since the potential has compact support, we can define the 
scattering length ao associated to V(x) by the equation f(x) = 1 — ao/|x| for x outside the support 
of V(x) (this definition can be generalized by ao = linv-i-oc r(l — f(r)), if V has unbounded support 
but still decays sufficiently fast at infinity). Another equivalent characterization of the scattering 
length is given by the formula 

J dxV(x)f(x) = 8vra . (3) 

Physically, ao is a measure of the effective range of the potential V{x). 
The Hamiltonian of the A-boson system is then given by 

N 

H = -J2&j + T,V a (x i -x j ) (4) 

j=l i<j 

with V a (x) = (ao/a) 2 y((ao/a)x). By scaling, V a has scattering length a. In the following we keep 
ao fixed (of order one) and we vary a with A, so that when A tends to infinity a approaches zero. 
In order for the Gross-Pitaevskii theory to be relevant we have to take a of order A -1 (see for 
a discussion of other possible scalings). In the following we choose therefore a = ao/N, and thus 
V a (x) = N 2 V{Nx). Note that, with this choice of a, the Hamiltonian (j3J can be viewed as a special 



2 



case of the mean-field Hamiltonian 



- 1 

H mi = - A,- + - f3 3 V((3( Xl - Xj )) . (5) 



i<j 



The Gross-Pitaevskii scaling is recovered when (3 = N . We study the dynamics generated by (J5J) for 
other choices of (3 (j3 = N a , with a < 3/5) in [3]. 

Since we have N particles in a box of volume one, the density is given by p = N. Hence, the 
total number of particles interacting at a given time with a fixed particle in the system is typically 
of order pa? ~ iV~ 2 <C 1 . This means that our system is actually a very dilute gas, scaled so that 
the total volume remains fixed to one. 

The dynamics of the TV-boson system is determined by the Schrodinger equation 

id t ipN,t = HipNj (6) 

for the wave function tpN,t £ L 2 (M. 3N , dx) . Instead of describing the quantum mechanical system 
through its wave-function we can describe it by the corresponding density matrix 7^ = {ipN^i^Njl 
which is the orthogonal projection onto ipNt- We choose the normalization so that TijNt = 1- The 
Schrodinger equation © takes the form 

idaN,t = [H, 7jv,t] • (7) 

For large N this equation becomes very difficult to solve, even numerically. Therefore, it is desirable 
to have an easier description of the dynamics of the system in the limit N — > 00, assuming we are 
only interested in its macroscopic behavior, resulting from averaging over the N particles. In order 
to investigate the macroscopic dynamics, we introduce the marginal distributions associated to the 

(k) 

density matrix jN,t- The A;-particle marginal distribution ^ N t is defined by taking the partial trace 

(k) 

over the last N — k variables. That is, the kernel of 7*^ is given by 



7^ ) t (x fc ;x' fc ) = / dx A r-fc7A^(xfc,x A r ifc ;x' fe ,x A r_ 



where 7jv,t(x;x') denotes the kernel of the density matrix Tjvt- Here and in the following we use the 
notation x = (xi, . . . ,xn), x fc = (^1, • • • , %k), x-N-k = (xk+i, ■ ■ ■ , aJjy), and analogously for the primed 
variables. By definition, the /c-particle marginal distributions satisfy the normalization condition 

Tr 7 }J] = l for all k = 1, . . . , N . 

In contrast to the density matrix jN,t, one can expect that, for fixed k, the marginal distribution 
7^ has a well defined limit 7^ for N — > 00 (with respect to some suitable weak topology), whose 
dynamics can be investigated. In particular, the Gross-Pitaevskii equation (^Q) is expected to describe 
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the time evolution of the limit 7^ of the one-particle marginal distribution, provided 7^ = \ut)(ut\ 
is a pure state. Eq. Q can be generalized, for 7^ describing a mixed state, to 

id t7 £] t (x;x') = (-A + A') 7 £ i (x;x') + Svrao (j£ >t (x;x) - ^{x'^x')) >y£] t (x; x') , (8) 

which we again denote as the Gross-Pitaevskii equation. 

To understand the origin of (JHJ), we start from the dynamics of the marginals 7!^+, for finite N. 
From the Schrodinger equation we can easily derive a hierarchy of N equations, commonly called 

(k) 

the BBGKY hierarchy, describing the evolution of the distributions "f N t , for k = 1, . . . , N: 



+ (N -k)^2 J &Ek+i(y a (xj - x k+1 ) - V a (xj - a;fc+i))7^ i t 1) ( x fe' a; fc+i; x fc' a; fc+i) 



i=i 

(k) (1) 

Here we use the convention that *y N t = 0, for k > N. Hence, the one-particle marginal density 7^ 
satisfies 

ida^ixr, x[) = {-A X1 + A a ./ i )7{J ) t (xi;xi) 

f ' (2) ( 10 ) 

+ (N-1) / dx 2 (Vaixx - x 2 ) -V a (x[ - x 2 )) T N > t (x 1 ,x 2 ;x' 1 ,x 2 ) ■ 

In order to get a closed equation for 7]^ we need to assume some relation between 7^ and 7^. 
The most natural assumption consists in taking the two particle marginal to be the product of two 

identical copies of the one particle marginal. Although this kind of factorization cannot be true for 

(k) 

finite N, it may hold in the limit N — ► 00. We suppose therefore that "y^ t, for k = 1,2, is a limit 

Ik) 

point of J N t , with respect to some weak topology, with the factorization property 

(2) / / / \ (l) / / \ (l) / / \ 

Too,t{x 1 ,x 2 ;x 1 ,x 2 ) = rj^{x 1 ;x 1 )^ t {x 2 ;x 2 ) . 

Under this assumption we could naively guess that, in the limit N — > 00, Eq. (|1U|) takes the form 

idaS^x',) = (— A X1 + A 4 ) 7 « + (Q t (x x ) - Q t &)) l^i^'i) ( n ) 

with 

Q t (xi)= lim N / dx 2 - x 2 )Y^ t (x 2 ; x 2 ) 

= lim / dx 2 N 3 V(N( Xl -x 2 )) 1 £ >t (x 2 ;x 2 ) 
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where we defined bo = J dxV(x). Using the last equation, can be rewritten as 

ida^xr,^) = (-A X1 + A^g^;^) + b (tS^i^i) -TS^'l^i)) (12) 

which is exactly the Gross-Pitaevskii equation (JHJ), but with the wrong coupling constant in front 
of the non-linear term (bo instead of 87rao). The fact that we get the wrong coupling constant 
suggests that something was not completely correct with the naif argument leading from 1)10(1 to 
Reconsidering the argument, the origin of the error is quite clear: when passing to the limit 

N — > oo we first replaced ^ N t with 7^ t and only after this replacement we took the limit N — > 00 

(2) 

in the potential. This procedure gives the wrong result because the marginal distribution 7]^ t 
has a short scale structure living on the scale 1/N, which is the same length scale characterizing 
the potential V a {x). The short scale structure of ^ Nt (which describes the correlations among the 

particles) disappears when the weak limit is taken, so that 7^^ lives on a length scale of order one. 

Therefore, in (|12|) we get the wrong coupling constant because we erroneously disregarded the effect 

(2) 

of the correlations present in 7jy- t . It is hence clear that in order to derive the Gross-Pitaevskii 

equation (JHJ) with the correct coupling constant 87rao, we need to take into account the short scale 
(2) 

structure of ^ N t - To this end we begin by studying the ground state of the system. 

A good approximation for the ground state wave function of the iV boson system is given by 

N 

W(x) = l[f(N(x i -x j )) 

i<j 

where f(x) is defined by (J2J) (then f(Nx) solves the same equation (J2J) with V replaced by V a ). Since 
we assumed the potential to be compactly supported (let R denote the radius of its support), we 
have f(x) = 1 — ao/|x|, for \x\ > R, and thus f(Nx) = 1 — clq/N\x\ = 1 — a/|x|, for |x| > Ra. 
A similar ansatz for the ground state wave function was already used by Dyson in [2] to prove his 
upper bound on the ground state energy. In order to describe states of the condensate, it seems 
appropriate to consider wave functions of the form 

Vw(x) = W(x)4> n (xl) 

where </>at(x) varies over distances of order one, and is approximately factorized, that is </>jv(x) ~ 
n^Li 0( a; j)- Assuming for the moment that this form is preserved under the time-evolution we have 

jP jt (xi,x 2 ;x' 1 ,x' 2 ) ~ f(N(xi -x 2 ))f(N(x' 1 - x' 2 ))^] t {xi] x[)^ t (x2; x' 2 ) . 
(2) 

Thus, for finite N, jj^ t is not exactly factorized and has a short scale structure given by the function 
f(Nx). When we consider the limit N — » 00 of the second term on the right hand side of (|1U() we 
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obtain 

„(2) , 



lim N / dx 2 V a (x 1 - x 2 )T N , t (x 1 ,X2;x' 1 ,x 2 ) 



lim N 3 [ dx 2 V(N( Xl - x 2 ))f(N( Xl - x 2 )h£] t (x 1] x[h£] t (x 2 ;x 2 ) ( 13 ) 



= 8™oTw,t( x r, x ihS,t( x r, x i) 

where we used Eq. (j3J) and the fact that 7SI lives on a scale of order one (and thus we can replace 
it by 7^ without worrying about the correlations). This leads to the Gross-Pitaevskii equation for 

/oo,t> 

^StOl^'l) = (- A x! + A^J 7^(xi;x / 1 ) + 8vra (j^ t (xi;xi) - 7^(3^ 5 x i)) lS,t( x ^ x l) 

which has the correct coupling constant in front of the non-linear term. 
Note that the factorization 

T^,t( x u x ^ x 'i, x 2) = lS,t( x ^ x ihS,t( x 2-,x' 2 ) 

(2) 

still holds true, because the short scale structure of j N t vanishes when the weak limit N — > 00 is 
taken. The short scale structure only shows up in the Gross-Pitaevskii equation due to the singularity 
of the potential. 

In order to make this heuristic argument for the derivation of the Gross-Pitaevskii equation 
rigorous, we are faced with two major steps. 

i) In the first step we have to prove that the fc-particle marginal density in the limit N — > 00 
really has the short scale structure we discussed above. That is we have to prove that, for large 
N, ' 

/fc+i \ 
7$f Wi;4+l) ~ f J] f(N( Xi - Xj ))f(N(x' t - x' 3 ))\ 7 2j 1) (x fc+ i; x' fc+1 ) (14) 

where 7^^ ^ s * ne limit of J^ 1 ^ with respect to some suitable weak topology (in the heuristic 
argument above we considered the case k = 1, here k is an arbitrary fixed integer k > 1). Eq. 
(f 141) would then imply that, as N — > 00, the last term on the r.h.s. of the BBGKY hierarchy 
@ converges to 

lim N / dx fc+ iK(xj -x fc+ i)7^ 1} (x fe+1 ;xfc +1 ) = 87ra 7^t 1) (x fe ,x i ;x / fc ,x j )- 

TV — »oo J ' 
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Therefore, if we could also prove that the second term on the r.h.s. of @ vanishes in the limit 
N — > oo (as expected, because formally of the order iV _1 ), then it would follow that the family 
satisfies the Gross-Pitaevskii hierarchy 

k 



3=1 

+ 87ra J dxfc+i (5(x k+1 - Xj) - 5(x k+1 - x'j)) 7^ 1) (x fe , x k+1 ; x' fc , 



3=1 

for all A; > 1. We already know that this infinite hierarchy of equation has a solution. In fact 
the factorized family of densities 7^ t (xfc; x' k ) = IJj=i 7co!*( x i' ^j) ^ s a solution of (JT5|) if and 
only if 7^ solves the Gross-Pitaevskii equation Q . 

ii) Secondly, we need to prove that the densities 7^ factorize, that is, that, for all k > 1, 

k 

r£t ( x fc ; 4 ) = II 7 S* ( x i ]X '^- ( 16 ) 

i=l 

Then, from l|15|) and (|16|) . it would follow that 7^ is a solution of the Gross-Pitaevskii equation 
©• Note that, since we already know that (|15[) has a factorized solution, in order to prove 
(|16j) it is enough to prove the uniqueness of the solution of the infinite hierarchy (|15[). 

Unfortunately, due to the singularity of the <5-function, we are still unable to prove that ()15|) has a 
unique solution and thus we cannot prove part ii) (the best result in this direction is the proof of the 
uniqueness for the hierarchy with a Coulomb singularity, see jS] ) . On the other hand we can complete 
part i) of our program, that is, we can prove that any limit point {7^ \}k>l °f the family {7^l}feLi 
(with respect to an appropriate weak topology), satisfies the infinite hierarchy ()15|) . provided we 
replace the original Hamiltonian H with a slightly modified version H, where we artificially modify 
the interaction when a large number of particles come into a region with diameter much smaller than 
the typical inter-particle distance. Since H agrees with H, apart in the very rare event (rare with 
respect to the expected typical distribution of the particles) that many particles come very close 
together, we don't expect this modification to change the macroscopic dynamics of the system: but 
unfortunately we cannot control this effect rigorously. 

Note that the Gross-Pitaevskii equation Q is a nonlinear Hartree equation 

id t ut = -Au t + (V * \u t \ 2 )u t (17) 

in the special case V(x) = 8Trao5(x). In the literature there are several works devoted to the 
derivation of (|17|) from the iV-body Schrodinger equation. The first results were obtained by Hepp 
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in [5], for a smooth potential V(x), and by Spohn in ^5], for bounded V(x). Later, Ginibre and 
Velo extended these results to singular potentials in 0: their result is limited to coherent initial 
states, for which the number of particles cannot be fixed. In [Sj, Erdos and Yau derived (|17[) for the 
Coulomb potential V(x) = ±l/|x|. More recently, Adami, Bardos, Golse and Teta obtained partial 
results for the potential V(x) = 5(x), which leads to the Gross-Pitaevskii equation, in the case of 
one-dimensional systems; see [Q. 



2 The Main Result 

In this section we explain how we need to modify the Hamiltonian and then we state our main 
theorem. In order to derive (|15|) it is very important to find a good approximation for the wave 
function of the ground state of the N boson system. We need an approximation which reproduces 
the correct short scale structure and, at the same time, does not become too singular (so that error 
terms can be controlled). Our first guess 

W(x) = J] f a ( Xi - Xj ) = n f(N( Xi - Xj )) (18) 

is unfortunately not good enough. First of all we need to cutoff the correlations at large distances 
(we want f a (x) = 1 for \x\ ^> a). To this end we fix a length scale i\ ^> a, and we consider the 
Neumann problem on the ball {x : \x\ < £\} (we will choose i\ = jV _2 / 3+K for a small k > 0). We 
are interested in the solution of the ground state problem 

(-A + l/2V a (x))(l - w(x)) = e h (1 - w(x)) 

on {x : \x\ < li}, with the normalization condition w(x) = for |x| = t\. Here is the lowest 
possible eigenvalue. It is easy to check that, up to contributions of lower order, ei x ~ ?>a/l\ . We can 
extend w(x) to be identically zero, for \x\ > i\. Then 

(-A + l/2V a (x))(l - w{x)) = q(x)(l - w(x)), with 

q( x ) - ^x(\x\ < h) ■ 19 ^ 

For a <ti \x\ <^ t\, the function 1 — w(x) still looks very much like 1 — a/|sc| 3 but now it equals one, 
for |x| > l\. Replacing f a {xi — Xj) by 1 — w(xi — Xj) in ((TH|) is still not sufficient for our purposes. 
The problem is that the wave function ni<j(l — w i x i ~ x j)) becomes very singular when a large 
number of particles come very close together. We introduce another cutoff to avoid this problem. 
We fix a new length scale I S> t\ S> a, such that I iV -1 / 3 (that is t is still much smaller than 
the typical inter-particle distance: we will choose i = iV~ 2 / 5 ~ K for a small k > 0). Then, for fixed 
indices i and j, and for an arbitrary fixed number K > 1, we cutoff the correlation between particles 
i and j (that is we replace 1 — w(xi — Xj) by one) whenever at least K other particles come inside 
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a ball of radius t around % and j. In order to keep our exposition as clear as possible we choose 
K = 1, that is we cutoff correlations if at least three particles come very close together. But there 
is nothing special about K = 1: what we really need to avoid are correlations among a macroscopic 
number of particles, all very close together. To implement our cutoff we introduce, for fixed indices 
i,j, a function -F^x) with the property that 

i^-(x) 1 if j \ Xi ~_ Xm \ ^ 1 for all m + i,j 

{ | Xj x m | ^> t 

Fij(x) = otherwise. 

Instead of using the wave function n«<j(l — w ( x i ~ x j)) we wm approximate the ground state of the 
N boson system by 

W(x) = Y[(l-w(x i -x j )F ij (x)) (20) 

i<j 

(the exact definition of W^(x) is a little bit more complicated; see @], Section 2.3). The introduction 
of the cutoffs Fy in the wave function W(x) forces us to modify the Hamiltonian H. To understand 
how H has to be modified, we compute its action on W(x). We have, using (|19|). 

iy(x)- 1 (W)(x) = Y^qixi-Xj) 

id 

+ Y {{l/2)V a ( Xi - Xj) - q(x t - Xj)) (1 - Fy(x)) 

+ lower order contributions. 

The "lower order contributions" are terms containing derivatives of Fij : they need some control, but 
they are not very dangerous for our analysis. On the other hand, the second term on the r.h.s. of the 
last equation, whose presence is due to the introduction of the cutoffs Fij, still contains the potential 
V a and unfortunately we cannot control it with our techniques. Therefore, we artificially remove it, 
defining a new Hamiltonian H, by 

H = H-Y J ((l/2)V a (xi - x 3 ) - q(xi - Xj)) (1 - i^(x)) . 

Note that the new Hamiltonian H equals the physical Hamiltonian H unless three or more particles 
come at distances less than i <^ A r ~ 1 / 3 . This is a rare event, and thus we don't expect the modification 
of the Hamiltonian H to change in a macroscopic relevant way the dynamics of the system. 

Before stating our main theorem, we still have to specify the topology we use in taking the limit 
N — > oo of the marginal distributions 7jyl- It is easy to check that, for every k > 1, 7^j(xfc;x / fc ) € 
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L 2 (A k x A fc ). This motivates the following definition. For T = {7^}fc>i G fc>1 L 2 (A fe x A fc ), and 
for a fixed v > 1, we define the two norms 

||r||_ := J>~ fc ||7 (fe) l|2 and ||r||+ := sup i/ fc || 7 (fc) || 2 (21) 
fc>0 fc - 1 

where ||.||2 denotes the L 2 -norm on A fc x A fc . We have to introduce the parameter v > 1 to make sure 
that, for Tj^ jt = {^\}^ = i, the norm [|rjv,t||- is finite (choosing v large enough, we find ||rjv,t||- < 1, 
uniformly in N and t). We also define the Banach spaces 

W- := {r = {7 ( * ) }*>0 £®i 2 (A* x A fe ) : ||r||_ < oo} 

fc>l 

and 

H + := {T = { 7 (/c) H>o G©i 2 (A fc x A fc ) : lim ^||7 (fc) || 2 = 0}. 

k>l 

Then we have (Ti.-, ||-||-) = (T~C+, ||.||+)*. This induces a weak* topology on with respect to 
which the unit ball £>_ of 7i- is compact (Banach- Alaouglu Theorem). Since the space 7Y+ is 
separable, the weak* topology on the unit ball B- is metrizable: we can find a metric p on TL-, such 
that a sequence T n G B- converges with respect to the weak* topology if and only if it converges 
with respect to the metric p. For a fixed time T, we will also consider the space C([0,T],B-) of 
functions of t G [0, T], with values in the unit ball i3_ C 7~t~, which are continuous with respect to 
the metric p (or equivalently with respect to the weak* topology of H-). We equip C([0, T], 
with the metric 

£(ri(*),r 2 (t)) = sup p(ri(*),r 2 (*)). 

te[o,T] 

In the following we will consider the families Tpf,t = {7jvi}ib=i as elements of C([0, T], £>_), and we 
will study their convergence and their limit points with respect to the metric p. We are now ready 
to state our main theorem. 

Theorem 2.1. Assume a = clq/N, l\ = N~ 2 ' 3+K , i = N~ 2 ^ 5 ~ K , for some sufficiently small k > 0. 
Assume 

(ij N ,o,H 2 ^ Nfi )<CN 2 , 

where (.,.) denotes the inner product on L 2 (M 3Ar ,dx). Let ipN,t, fort G [0, T], be the solution of the 
Schrodinger equation 

id t ^ N ,t = H1>N,t (22) 

with initial data ipN,o- Then, if a = (||V||i + || V||oo) is small enough (of order one) and v > \ is 
large enough (recall that v enters the definition of the norms \21))). we have: 

i) ^N,t = {7jvj}fcLi has at least one (non-trivial) limit point = { 7 ^\}fc>i G C([0,T],B-) 
with respect to the metric p. 
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ii) For any limit point = {7^\}fc>i and for all k > 1, there exists a constant C such that 

Tr(l-A l )(l-A J ) 7 S<C (23) 

for all iytj,t£ [0,T]. 

in) Any limit point satisfies the infinite Gross-Pitaevskii hierarchy h!5\) when tested against 
a regular function j( fc )(xfc; x' fc ): 

-8«7ra o y] / ds / dx fc dx' fc J (fc) (x fc ; x' fe ) ( 24 ) 
x / dx fc+ i (5(xj - x k+ i) - 5(x'j - x k+1 )) 7^+ 1) (x fc ,x fc+ i;x / fc ,x A:+ i) . 



Here we use the notation (J^ k \-y^) = J dx^dx'^ J( fc )(x^; xC) ^^(x^; x' fc ). 

i) The main assumption of the theorem is the requirement that the expectation of H 2 at t = 
is of order N 2 . One can prove that this condition is satisfied for ^jv,o( x ) = W(x)c/>n ( x ) and 
i^at sufficiently smooth (see [3], Lemma Dl). Physically, this assumption guarantees that the 
initial wave function V>tv,o( x ) has the short-scale structure characteristic of W(x) and, hence, 
that it describes, locally, a condensate. 

ii) It is a priori not clear that the action of the delta-functions in the Gross-Pitaevskii hierarchy 
(|24[) is well defined. This fact follows by the bound (|23j) . which makes sure that is 
sufficiently smooth. 

iii) We also need to assume that a = ([|V|| o0 + ||V||i)is small enough (but still of order one). This 
technical assumption is needed in the proof of the energy estimate, Proposition 13.11 

3 Sketch of the Proof 

In this section we explain some of the main ideas used in the proof of Theorem 12.11 Let ipN,t be the 
solution of the Schrodinger equation ()22|) (with the modified Hamiltonian H). We can decompose 

ipN,t as 
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where W(x) is the approximation for the ground state wave function defined in (|2U|), This decom- 
position is always possible because W(x) is strictly positive. 

The main tool in the proof of Theorem l2.1l is an estimate for the L 2 -norm of the second derivatives 
of (j>N,t- This bound follows from the following energy estimate. 

Proposition 3.1. Assume a = clq/N, l\ = _/V~ 2 / 3+K and I = N~ 2 / 5 ~ K for k > small enough, and 
suppose a = (||V||i + H^Hoo) is sufficiently small. Then there exists a constant C > such that 

N 

dx|(W0)(x)| 2 > (C-o(l)) Y, / dxTy 2 (x)|V,V^(x)| 2 

-o(l) (^NJ2 J dx^ 2 (x)|V l 0(x)| 2 + iV 2 y dxVF 2 (x)|0(x)| 2 ^ , 

where o(l) — > as N — > oo. 

Remark. The proof of this proposition is the main technical difficulty in our analysis. It is in 
order to prove this proposition that we need to introduce the cutoffs in the approximate ground 
state wave function W(x), and that we need to modify the Hamiltonian. 

Using the assumption that, at t = 0, (ipN,o,H ipNfi) < CN 2 , the conservation of the energy, and 
the symmetry with respect to permutations, we immediately get the following corollary. 

Corollary 3.2. Suppose the assumptions of Provosition W^7\ are satisfied. Suppose moreover that the 
initial data tpN,o is symmetric with respect to permutations and (?/>jv,o> -f^V'N.o) <■ CN 2 . Then there 
exists a constant C such that 

y W 2 (x)|V i V J ^, i (x)| 2 <C (25) 
for all i 7^ j, t and all N large enough. 

Remark. The bound (|25|) is not an estimate for the derivatives of the whole wave function ipN,t- 
The inequality 

J dx|ViV^(x)| 2 < C (26) 

is wrong, if ip satisfies (ip,H 2 ip) < CN 2 . In fact, in order for (■ip,H 2 ip) to be of order TV 2 , the 
wave function V( x ) needs to have the short scale structure characterizing W(x). This makes ()26(l 
impossible to hold true uniformly in N. Only after the singular part W(x) has been factorized out, 
we can prove bounds like (|25f) for the derivatives of the remainder. One of the consequences of our 
energy estimate, and one of the possible interpretation of our result, is that the separation between 
the singular part of the wave function (living on the scale 1/N) and its regular part is preserved by 
the time evolution. 
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Next we show how the important bound ()25j) can be used to prove Theorem 12.11 According to 
the decomposition ^>jv,t( x ) = W(x)</>jv,t(x), we define, for k = 1, . . . ,N, the densities J7^(x^; x' fc ), 
for = 1, . . . , N, to be, roughly, the /c-particle marginal density corresponding to the wave function 
4>N,t (the exact definition is a little bit more involved, see [I], Section 4). The estimate ()25|) for the 

(k) 

second derivatives of 0jv,t translates into a bound for the densities U N t : 

TtCl-AOCl-A^tfJJO <C (27) 

for alii, j < TV with i ^ j, for all i and for all JV large enough. 

Moreover, we can show that, for v > 1 large enough (recall that the parameter v enters the 
definition of the norms l|21jl). the families Ujyj = {U^\}^ =1 define an equicontinuous sequence in the 
space C([0,T],B-) (this follows from a careful analysis of the BBGKY hierarchy associated to the 
Schrodinger equation (|22|): see [H, Sections 9.1 and 9.2 for more details). Applying standard results 
(Arzela-Ascoli Theorem), it follows that the sequence £/jv,t has at least one limit point, denoted 
Uoo,t = {uf£\}k>i, m the space C([0, T], B-). The bound (f2*7|) can then be passed to the limit 
N — > oo, and we obtain 

Tr(l-A l )(l-A J )C/S, ) t <C 

for all i ^ j and i G [0,T]. 

Next we go back to the family rV,t = {7jv+}fc=i- Clearly, the densities do not satisfy the 
estimate (|27|), In fact, still contains the short scale structure of W(x) (which, on the contrary, 

has been factorized out from U^\), and thus cannot have the smoothness required by (|27jl. 

It is nevertheless clear that the short scale structure of Tpf,t disappears when we consider the 
limit — > oo (in the weak sense specified by Theorem 12. 1|) . In fact, one can prove the convergence 
of an appropriate subsequence of T^t to the limit point i/oo.t of J7jv,t- In other words one can show 
that limit points of T^ jt , denoted by Tooj, coincide with the limit points of U^,t- Therefore, even 
though Tpfj) fo r finite N, does not satisfies the bound (|27|). its limit points = {j^\}k>i do. 
For every k > 1 we have 

Tr(l-A 4 )(l-A J )7S<C (28) 

for all i ^ j and t S [0, T]. This proves part i) and ii) of Theorem 12 . II fthe non-triviality of the limit 
follows by showing that Tr 7^ = 1). The bound (|28|) can then be used to prove part iii) of Theorem 
12.11 that is to prove that the family satisfies the infinite Gross-Pitaevskii hierarchy (|24jl. In 

(k) 

fact, having control over the derivatives of 7^, t allows us to prove the convergence of the potential 
to a delta-function (that is, it allows us to make (|13|) rigorous). To this end we use the following 
lemma (see jl], Section 8). 

Lemma 3.3. Suppose 6p(x) = f3~ 3 h(x/f3), for some regular function h, with f h{x) = 1. Then, for 
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any 1 < j < k, and for any regular function J(x k ;x' k ), we have 

j dxfcdx^dxfc+i J(x k ;x' k )(5p(xj - x k+1 ) - 5(xj - x k+ i))^ k+1) (x fe , x k+1 ; x k , x k+1 ) 

<Cj^P Tr(l - A,)(l - A fe+1 ) 7 ( fe+1 ) . 

Part iii) of Theorem 12.11 can then be proven combining this lemma with the estimates (|27|) and 
((2*5|) (see [I], Section 9.4, for more details). 
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